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Abstract 

An integro-differential equation satisfied by an eigenvalue density denned as the loga- 
rithmic derivative of the average inverse characteristic polynomial of a Wilson loop in 
two dimensional pure Yang Mills theory with gauge group SU(N) is derived from two 
associated complex Burgers' equations, with viscosity given by 1/(2N). The Wilson loop 
does not intersect itself and Euclidean space-time is assumed flat and infinite. This 
result provides an extension of the infinite N solution of Durhuus and Olesen to finite 
N, but this extension is not unique. 
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1 Introduction. 

In [1], in the context of two dimensional YM theory on the infinite plane with gauge 
group SU(N), the function 4>n(,U,t) defined by 

<Pn(v, r) = log [e^ii-y) (det(e» + W))] (1) 



was shown to satisfy Burgers' equation 



1 d 2 - 



with initial condition 

<j>N{y,o) = _ 2 tanh | ^ 3 ' ) 

is a Wilson operator associated with a non-selfintersecting loop, r measures the area 
enclosed by the loop in units of the 't Hooft gauge coupling. (...) denotes averaging with 
respect to the exponent of the two dimensional YM action, r > and y and <^jy are 
real. 

At N = oo a shock appears at r = 4; for finite iV the shock is smoothed in a universal 
way in the regime y ~ 0, r ~ 4. (f>N admits a pole expansion with exactly integrable pole 
dynamics in a "time" r. 

The N = oo critical value r = 4 corresponds to the Durhuus-Olesen (DO) [2j phase 
transition point. That phase transition was found by solving a complex inviscid (N = oo) 
Burgers' equation. More specifically, y was taken to approach purely imaginary values 
and 4>n was complex, r remained real and non-negative. 
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The objective of this letter is to identify (det(z — W) -1 ) as an object that is more 
directly linked to the DO solution. In principle, there exist an arbitrary number of 
possible extensions of the DO equation and solution; this choice is special because all 
finite N effects are accounted for by a viscous term in an associated Burgers' equation 
with viscosity given by similarly to pp. The singularity structure for finite N is 
richer here than in pp. 



2 Conventions. 

The partition function of Euclidean 2D SU (N) YM is written as: 



, ——j — f SxtrF^F^ 

Z= [VAJe ^ J (4) 



't Hooft's coupling is 

A = g 2 YM N (5) 

and the area enclosed by the loop is A. 

The integration over at fixed W induces a probability density for W, given by 

V N (W,XA) = J2 d RXR(W)e-M c *W ( 6 ) 

R 

where R denotes an irreducible representation of SU (N) of dimension cLr and the char- 
acter Xr(W) satisfies: 

XR (1) = d R , J dW X R(W)x*s(W) = 8 RS (7) 

dW is the normalized Haar measure on SU(N) and Vn is defined relative to it. C2(R) 
is the quadratic Casimir of R; for the defining representation F we have: 

C 2 (F) = N-± (8) 

Vn obeys: 

V N (W, XA) = V N (W*, XA) = V* N (W, XA) = V N (W\ XA) (9) 
The variable r above is given by 

r = XA (l + 1) (10) 

The probability density can be viewed as a function of r: 

Pn(W,t)=V n (W,XA) (11) 
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3 Antisymmetric representations. 
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Expanding the characteristic polynomial we obtain a sum over all /c-fold antisymmetric 
representations F Ak of dimension dk = (^): 



det(z - W) = z N 
For W = diag(e l611 , e l6,JV ) we have 



Also, 



Xk (W) = jtoi +e * + "" +e iJ 
l<3i<32--<3k<N 



C 2 (k) = ^p-HN -k) = C 2 {N- k) 



(12) 



(13) 



(14) 



For the computation of (det(z-W)) only the representations R = F Ak in the sum giving 
Pjv contribute. One obtains, with real y: 

<det(e» + W)} = e f(^"i)^ ^ e f( fc -f ) e 2w( fe -f) 2 = e f (^V(y) (15) 

fc=o V fc / 

The forms of the prefactor of y in the exponent and the prefactor of r in the exponent 
show that qN{y) obeys the (linear) heat equation 



dqN 1 d 2 q 



N 



dr 2N dy 2 



(16) 



which leads to ([2]) 



4 Symmetric representations. 



Equation (I16p holds as a consequence of the linearity of P/v in representation space, the 
restriction of contributions to a subset that can be labeled by an index k, and the 
quadratic dependence of C2(Rk) on k. The dimensions d^ only enter through the initial 
condition. Therefore, one expects a similar equation to hold for a sum over all fc-fold 
symmetric representations. This time the range of k extends to infinity, the generating 
function of all the R^ is no longer a polynomial and, consequently, a richer analytic 
structure is expected. 

The generating function for all symmetric representations is the inverse of the char- 
acteristic polynomial. The general formula can be obtained by considering a diagonal 
W: 
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det(z - W) 



1 

7n 



N V N 



n 



j=l ln=0 



E 



z -n e in9j 



1 

7N 



k=l 



Xk{W) 



(17) 
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where the character and dimension of each Rk are given by: 



Xk (W)= e ime 1 +in 2 02+....+in k h . dk 

?ii,n2,...,njv>0, y^'^_ 1 rij=k 



N + k-r 
N-l . 



(18) 



Most importantly, the second order Casimirs are quadratic in k: [3] 

C 2 (k) = ^-±k(N + k) 



N 



A new area variable, t, now replaces r: 

N-l 



t = t- = XA[1 

N + l V N 



(19) 



(20) 



Qn{u,t) is replaced by two functions of z, ip^. y> (z, t). i s defined for \z\ > 1 and 

is defined for \z\ < 1. The ip± are analytic in their respective domains. 



(det(2- W)- 1 ) =ip { ±\z,t) 



(21) 



+ or — hold, depending on whether z is inside the unit circle or outside it. For \z\ > 1, 



ip±^\z,i) are given by: 



z s z 



At t = Pj\r(W,0) = S(W, 1) with respect to the Haar measure. Therefore, 

1 



#°M) 



(z-l 



The linear equation replacing (|16p is: 



0f 



2N\ dz 2 ^ ± v ' 7 ' 



depending on the domain of z. <j)( N \y,T) is replaced by: 



4*° cm) 



5 iV 



Explicitly, 



1 1 Eg°=i fc4^ Tfc e- f ^- 

"T" »t k(N+k) 

2JV 



2 ^l + EfeLi dfc^e-* 



These functions obey 



^4)>(M) + ( i2 - 



4 W (*.«) 1 

<9z 



at 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



As before, an exponential substitution leads to Burgers' equation: 

z = e' lY (28) 

The map Y — > z takes the real axis into the unit circle, the ^sY > half plane into 
\z\ > 1 and the QY < half plane into \z\ < 1. Every strip \?RY — 2kir\ < it, k G Z 
is mapped onto the entire z-plane. In terms of Y, there is only interest in functions 
periodic under Y — > Y + 2kn, which define single valued functions of z. Viewing the 
as functions of Y the nonlinear PDE-s become complex Burgers' equations, 

2N~dY^ = ~dT + ~dY~^ (29) 

#)( e -^,0) = icot| (30) 

The explicit forms of the solutions in terms of sums over k (|22l26p imply the following 
asymptotic behavior at t — > oo: 

^\e- iY ,oo) = 4 (31) 



with initial conditions: 



5 Relation to DO. 

At N = oo, ^°\e- iY ,t) is related to the function f(A, a) of [2] by: 

f{A,a) = ^+ ) (e- iy ,t) (32) 

with 

Y = a, t = — (33) 

Therefore, ^ {e~ iY , t) is one possible extension of the DO solution to finite N. 
At infinite N DO define the eigenvalue density of W, pa(cx), now for real a, by 

p A (a) = -2 lim [3/(A,a)] (34) 

Sa— +0+ 

Observe that for any i > the definitions of (z, t) by sums over k can be analytically 

extended to all z, excepting z = for ip^P and z = oo for ip_. In particular, for t > 0, 

(z, t) are well defined for |z| = 1. The unit circle \z\ = 1 is parametrized by z = e~ ty 
with real y. 

It is easy to check that 0+ (e~ iy , t) + <j)_ {e~ %y ,t) is purely real and 4> ( ^\e~ w , t) — 

c()^ N \e~' iy , t) is purely imaginary. This finally leads to an extension of the DO infinite N 
eigenvalue density to finite N: 

P {N) (y,t) = ±[<f>W(e-*,t) t)\ (35) 
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The limiting behavior at t = oo is now seen to be 

1 

2^' 



(y,oo) 



(36) 



and is N independent. The N independent initial condition also requires a limiting 
procedure because singularities appear at z = 1 when t attains the value 0: 



P (7V) (y,0) = ^lim t ^ 0+ {lim £ 



+0+ 



(AO, 



-ly-e 



,*)]} 



2tt 



lirm 



cot 



y—it 



cot 



,<5(j/-2fc7r) 



(37) 



The initial condition is also iV-independent. Thus, the entire iV-dependence of pW is 
contained in the differential equations, more specifically, in their viscous terms. 



6 Equation for p( N \ 



At infinite N p^ is the Wilson loop matrix eigenvalue density and therefore a more 
physical object than the average of the inverse characteristic polynomial. It therefore 
seems desirable to derive an equation for p( N ) (y, t) directly, without the involvement of 
other functions. 

The equations obeyed by 4^± \e~~' iy , t) have only one nonlinear term 



ia_ 

Td~y 



which hinders superposition. Using 



, f .y 



one could get an equation just for u 



CAT) _ W 
IN) AN) 



(N) + AN) 



(38) 



(39) 







(N) + JN) 



(e %y ,t) in terms of the difference 



(e iy ,t) if one expressed the sum 



(e iy , t). This is pos- 



.(N) AN) 



sible since, for t > 

k(N), 



,W _i_ AN) 



(e iy ,i) are the real and imaginary parts of the 



analytic function (f>\ (z,t) on the curve \z\ = 1. 

The needed device is the Hilbert transform H, mapping functions of a real variable, 
g{y) into other functions of a real variable, (Hg)(y): 



(Ug)(y) = lim - / 

A^oc 7T J- 



9{x 



■dx 



TT 



9{x) 



dx 



-AV -X 7T J__oo y - X 

H is defined with the opposite sign relative to [5]. For real y and t > one obtains: 

,2 / iaj\ . .... . \ 2 



(40) 



(^(e^tlj-PteAt))' 
The sought after equation follows: 



i(27r) 2 p( Ar )(y,t)(Hp^0(y,i) 



2N 



dy 2 



Of 



7T 



dy 



(y,t)(HpW)(y,t) 



(41) 



(42) 



This integro-differential equation is the main result of this letter. The equation is a 
particular case of an equation studied in [4]. The equation has been further investigated 
in [5]. The results of [5] were applied in the context of infinite N to 2D SU(N) YM 
by Blaizot and Nowak in [6j. Note that the nonlocal term already contributes in the 
inviscid limit. Thus, the result here provides a direct and minimal extension to finite N. 

Perhaps the most interesting property of this integro-differential equation is that 
turning on the viscosity does not assure regularity even for smooth initial conditions: 
one can have finite time "blow-ups", even for finite N. In our application we know 
that we start from a singular initial condition. Another interesting property of this 
equation is that it admits solutions given by superpositions of pole terms with the entire 
t dependence given by the location of the poles in the complex plane. The motion of the 
poles is governed by coupled first order differential equations of Calogero type. 

More work on the consequences of the above for physics is left for the future. 



7 Integral representations. 



For t > and \z\ ^ 1, equations (I22p admit an integral representation which includes the 
t —* + initial condition. The basic step is to write a Gaussian integral representation 
for the i-dependent term: 

r k(N+k) _ p ^ r dx i g 2 +i /T^+ff) 



e -w fc i JV+fc ) = e T / e ~2 x + ix Vn^+-> (43) 

D V27T 



So long as t > and \z\ > 1 the sum giving ip+^ can be interchanged with the integral 
and after that the sum over k can be performed. Changing variables of integration to 
u = x^J~jj produces: 

V>fW) =^^| f due-*" 2 (ze-^-e^ (44) 

From this, one can immediately get an integral representation for if)^ . For comparison, 
the analogue equation in the anti-symmetric case [I], for any z, is 



{det(z- W)) = e~ — W — J due~^ u [ze~ - el J (45) 

Eq. (|45p can be obtained from eq. (|44p with the formal replacements u — * —iu and 
y/N — > iy/~N, the latter causing also the replacement t — > r because we view A.A as fixed. 

These integral representations produce asymptotic expansions in 4= starting from a 
dominating saddle point and make evident the difference in analytic structures in z at 
finite N. 
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8 Discussion. 



The first choice for an analytic function containing information about the eigenvalue 
density of W would be the average resolvent: 

(R N (z,t)) = ±(tr(z-W)- 1 ) (46) 

i?jv has a natural expansion in n-wound loops, trW n , which can be converted to an ex- 
pansion in representations, but at the expense of linearity; all representations contribute 
to (R N ). 

Here and in pQ, I focused on the characteristic polynomial because: 

R N (z,t) = ±1£ log[det(z - W) ±l ] (47) 
TV oz 

For simple, properly normalized gauge invariant observables Oi one has infinite TV" fac- 
torization: 

<n^>=n<^> ( 4§ ) 

i i 

This implies that at infinite N 

(f((D)) = f((0}) (49) 

for an arbitrary function / at points where it can be Taylor expanded. Hence, (R^°°^ (z, t)) 
will be obtained with either sign choice in (|47p . 

Choosing + in (|47D and averaging with respect to Pjy produces a sum over N simple 
poles. However, the averaged eigenvalue density at finite N should be reflected by a cut 
running round the unit circle, completely segregating the interior of the unit circle from 
its exterior. This indicates that the choice of — in (|47p is a more appropriate extension 
of the infinite N eigenvalue density. 

This discussion makes it evident why the extension to finite N is non-unique. The 
preferred extension would perhaps be {RN(z,t)), but the method used for the charac- 
teristic polynomials in this letter fails there, and there is doubt that a simple finite N 
equation exists. There are other methods that yield the same result for the characteristic 
polynomial, and these methods might extend to Rn, starting from 

(detf^S) (50) 



Z2-W , 

and taking z\ — ► Z2 subsequently. I hope to explore this more in the future. 

The longer view is to recall that only the 2D problem can reduce to simple equations, 
while the main interest is to focus on the large N non-analyticity and its universal 
smoothing out at finite N, features which do appear to extend to higher dimensions 
El H2 \TU[ lllj . There, the main new ingredient is the need to renormalize. I think that 
the most convenient object to renormalize is the average characteristic polynomial, the 
topic of PQ , but would not rule out the average of the inverse characteristic polynomial as 
deserving more study in this context either. I hope to be able to provide a renormalized 
framework for dealing with the average characteristic polynomial of a Wilson loop in 3D 
and 4D some time in the future. 
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